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Abstract 



We study kicked quantum systems by using the squeezed state approach. 
Taking the kicked quantum harmonic oscillator as an example, we demon- 
strate that chaos in an underlying classical system can be enhanced as well as 
suppressed by quantum fluctuations. Three different energy diffusions are ob- 
served in the kicked quantum harmonic oscillator, namely, localization, linear 
diffusion, and quadratic diffusion. 
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I. INTRODUCTION 



A squeezed state is a generalized coherent state, which has a wide application in many 
branches of physics such as quantum optics, high energy physics, etc. Recent years have 
witnessed a growing application of the squeezed state to the field of nonlinear dynamics 
and chaos |I|-|12|. Since the squeezed state approach is a kind of approximation of quan- 



tum mechanics, it is called a semiquantum approach by people in this community. The 
time evolution of expectation values and fluctuations of the squeezed state is thus named 
semiquantum dynamics, or squeezed state dynamics. 

The main purpose of the squeezed state approach is to study how quantum fluctuations 
manifest themselves on classical trajectories. This approach starts directly from quantum 
systems with no reference to classical limit. In fact, it has been shown || that the squeezed 
state dynamics exits even for systems without a well-defined classical dynamics. Generally, 
the squeezed state approach simplifies the quantum version, and provides a complementary 
way to the semiclassical method. In some special cases (unfortunately, we still do not know 
under what kind of conditions), it gives us better results than the semiclassical method 
|T|-§],|7||| . Therefore, in addition to the semiclassical approach, the squeezed state approach 
is also a very useful tool to study the problem of classical-quantum correspondence. 

The squeezed state approach has proven to be very successful in studying dynamical 
systems ranging from integrable to many-body nonintegrable systems. Among many others, 
we name just a few examples here. In calculating the ground state energy for quantum sys- 
tem with potential V(q) = —Vq/ cosh 2 aq, Tsui discovered that the ground state energy 
obtained by the squeezed state approach is much closer to the exact ground state energy 
than that obtained from the WKB method. More recently, Pattanayak and Schieve [11] ap- 
plied this approach to a classically chaotic system, for which the WKB method completely 
fails. They have successfully calculated low-lying eigenenergies which agree within a few 
percent with the pure quantum (numerical) results. In addition to the ground state and/or 
lower excited state energies, the squeezed state approach can provide a way to obtain correct 
eigenfunctions. For instance, in applying the squeezed state approach to a harmonic oscil- 
lator, we can obtain not only the eigenenergies but also the eigenfunctions (including the 
ground state) exactly, as we shall see in Appendix A, whereas the semiclassical method fails 
to yield the exact eigenfunctions. As a further example of many-body systems, we would 
like to mention that we have used the squeezed state approach to the one-dimensional quan- 



tum Frenkel-Kontoroval model ||13|| , which is a nontrivial many-body, non-integrable system. 
Our results show that the squeezed state approach very nicely captures the feature of the 
quantum effect, namely, the standard map which determining the coordinates of the classi- 
cal ground state, is renormalized to an effective sawtooth map in the quantum case. The 
squeezed state results agree well with that of the quantum Monte Carlo method 

In this paper, we would like to apply the squeezed state to study generic behaviors of 
kicked quantum systems. It is well known that, in the development of quantum chaos, kicked 
quantum systems play a very important role. Prototypes of these kicked systems are the 
kicked rotator and kicked harmonic oscillator (KHO). They represent two different classes 
of dynamical systems. On the one hand, the kicked rotator obeys the Kolmogorov-Arnold- 
Moser (KAM) theorem. Classically, as the kick strength increases, invariant curves gradually 
break up. When the kick strength exceeds a critical value of K c = 0.9716..., the last invariant 
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curve disappears, and bounded chaos turns into global chaos, characterized by unbounded 
diffusion in the momentum direction |14]]. Quantum mechanically, the diffusion follows the 
classical one only up to a certain time, after which it is completely suppressed, thus leading 
to dynamical localization ||15|| . This phenomenon was connected to the Anderson localization 
P~6f1 , and was confirmed experimentally |Hj| . 

On the other hand, since a harmonic oscillator is a degenerate system, the KHO model is 
out of the framework of the KAM theorem, that is, diffusion can occur along stochastic webs 
for any small kick strength. As a matter of fact, the KHO model is not a toy model, it stems 
from a real physical system. It describes a charged particle moving in a magnetic field, and 
under the disturbance of a wave packet HIBfl . It has l| degrees of freedom. Classically, it 
depends on the ratio between frequency of the harmonic oscillator and that of the external 
kicks (see Eq.([34])), this system displays abundant structures in phase space such as crystal, 
quasicrystal, and stochastic webs [p^JOl . Moreover, since the phase space of the KHO is 
unbounded and cannot be reduced to a cylinder as in the case of the kicked rotator model, 
numerical investigation of the quantum KHO is much more difficult than that of the kicked 
rotator model. Therefore, contrary to the kicked rotator, only a few attempts have been 
made in the quantum KHO [p0| -|24|j. A general picture about the quantum behavior of the 
KHO is still lacking. 

In this paper, using the squeezed state approach, we are able to obtain the diffusion 
behavior of the quantum KHO model not only numerically but also analytically. The paper 
is organized as follows. In Sec. II we give a brief introduction of the squeezed state approach 
for purpose of self-containment In Sec. Ill, by using the KHO model, we shall demonstrate 
that quantum fluctuations will enhance chaos at a small perturbation regime, whereas it 
will suppress chaotic diffusion at a large perturbation regime. In Sec. VI, we shall study the 
diffusion and localization phenomena. In Sec. V, comparisons between the squeezed state 
results and the very few available quantum results will be given. We shall conclude our 
paper with discussions and remarks in Sec. VI. In Appendix A, we apply this approach to 
quantize the harmonic oscillator which exactly yields eigenenergies and eigenfunctions. In 
Appendix B, we outline our procedure of the pure quantum computation for KHO in some 
special cases. 



II. SQUEEZED STATE APPROACH 

The squeezed state approach starts from the time-dependent variational principle 
(TDVP) formulation 

5jdt{*(t)\ih--H\*(t)) = 0. (1) 

Variation with respect to (^(t)\ and \^(t)) gives rise to the Schrodinger equation, and its 
complex conjugate (see, e.g. Ref. [p5fl), respectively. The true solution may be approximated 
by restricting the choice of states to a subspace of the full Hilbert space, and finding the 
path along which the above equation is satisfied within this subspace. In the squeezed state 
approach, the squeezed coherent state is chosen as |^/(t)). In this manner, as we shall see 
below, in addition to the dynamics of centroid of wave packet, we will also have equations of 



3 



motion for the fluctuations, i.e., the spread of wave packet. Therefore, this approach enables 
us to study the effects of the quantum fluctuations on dynamical behavior. 

The squeezed state is defined by the ordinary harmonic oscillator displacement operator 
T>(a) acting on a squeezed vacuum state S(/3)\0): 

\aP) = V(a)S(fi)\0) 

Dipt) = exp (aa + — a*a) , (2) 
S{/3) = exp 



i(/% +2 -/Ta 2 ) 



a + and a are boson creation and annihilation operators which satisfy the canonical commu- 
tation relation [a, a + ] = 1. The coherent state is just operator T> acting on the vacuum state 
|0>, 

|a)=P(a)|0). (3) 

In terms of the number eigenstate \n), the coherent state can be written as 

/ 1 \ 00 a n 
|a) = exp --|a| 2 ) -=|n). (4) 

a is eigenvalue of creation operator, i.e., 

a\a) = a\a). (5) 

Denoting, 

a = \a\e iv } (6) 

and taking the integral over angle rj from to 2n on both sides of Eq. (^), we obtain the 
number eigenstate expressed in terms of the coherent state, 

|n) = — exp ( -\a\ 2 ) \a\- n Vri. di]e~ mr] \a) . (7) 
2n \2 J Jo 

The coordinate and momentum operators are defined as 

Q=^(a + + a). (8) 
Thus we have expectation values and variances 

q=(n)\mt)\) = Jl («•+«)■ (9) 
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Aq 2 = (*(t)\(q-q) 2 \*(t)\) = hG, 

a p * = (*(t)\(p- P ym)\) = n (^ + 4n 2 c) . (io) 



The canonical coordinates (G, IT) were introduced by Jakiw and Kerman [26] for the quantum 
fluctuations, and its relation with (3 in Eq. (^) is 



G = -\ cosh |/5| + — sinh 

u _i /3*-P sinh cosh |/?| 

2 \/3\ |cosh|/3| + ||Sinh|/3|| 2 ' 1 ' 

In the framwork of the squeezed state, the Heisenberg uncertainty relation becomes 

AqAp = |V1 + 16G 2 n 2 > ~. (12) 
In fact, the squeezed state \aj3) is equivalent to the Gaussian-type state 



W)> = J^yJI eX P {ftiPQ ~ W)) eX P ) |0> 

= e-*+\ap) i (13) 

where 

q = 1 - -1 + 2zn, 

e- 24 * = (cosh |/3| + A sinh |/3|) . (14) 

From the TDVP, we obtain dynamical equations for expectation values and quantum fluc- 
tuations, 

q = — • = -— (15) 
<9p ' dq ' 



^G = — , WI = -— , (16) 

which are canonical equations of motion. The dot denotes the time derivative and the 
Hamiltonian function H is given by definition 

H = m)\Hmt)). (17) 

These equations (0), and fllBD give us a simple and clear picture about the motion of the 
expectation values as well as the evolution of the quantum fluctuations, which are responsible 
for quantum diffusion. If the Hamiltonian consists of separate kinetic and potential terms 
such as 
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H=\p 2 + V{q), (18) 



then the Hamiltonian function can be written as 

2 



H = -p 2 + V(q) + hi K — + 2GU 2 ) + <J exp 



HI) 



l\V(q). (19) 



In the limit of h — 0, this Hamiltonian function reduces to the classical Hamiltonian. 

Initial conditions 

In order to solve the equations of motion ([15]) and (0), appropriate initial condition 
for variables (q,p) and (G, U) should be posed. In principle, the initial condition must be 
physically meaningful. Thus the following two conditions are generally selected. 

(1) Minimum uncertainty. The initial state |\l/(to)) should satisfy the condition of min- 
imum uncertainty. Because Go is always larger than zero, from the uncertainty principle 
[Eq. flTJ)], we have 

n(t ) = 0. (20) 

(2) Least quantum effect. We need to determine the initial value of G. This can be 
achieved by requesting minimization of H with respect to G, i.e. 

For instance, for a harmonic oscillator (Appendix A) we have G(t ) = 1/2lj . 



III. ENHANCEMENT AND SUPPRESSION OF CHAOS 

In this section, we would like to discuss the effect of quantum fluctuations on classical 
chaos. It is commonly argued that quantum fluctuations suppress classical chaos due to 
interference, while enhancing chaos due to tunneling. This is a rather qualitative argument. 
With the help of the squeezed state approach, we are able to do a quantitative analysis of 
quantum fluctuations. 

Up to now, there have been only a few works on this issue. In some models, suppression 
of chaos comes about, while in others enhancement takes place. Zhang and co-workers 
observed suppression of chaos in kicked spin and the kicked rotator systems Using a 

one-dimensional problem with a Duffing potential without any external perturbation, so that 
both classical and quantum behaviors are regular, Pattanayak and Schieve || demonstrated 
that the squeezed state behavior is chaotic, and concluded that the quantum fluctuations 
induce chaos. The effect of enhancement was also confirmed in a kicked double well model 
|| . We shall see later that the KHO model provides a prototype for studying these two 
effects. The enhancement and suppression can be observed by increasing the strength of the 
kicks. 

The Hamiltonian of the quantum KHO model can be written as [ITS] 
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H 



J + ^ 2 + V(q)&r, 



where 



$t= E 5{t-nT). 



(22) 



(23) 



Using the squeezed state as a trial wave function of Hamiltonian ([22]), one can readily obtain 



H = H + Hf + H, 



v 



where 



H 



P 2 , u 2 q 2 



2 2 

is the Hamiltonian of the harmonic oscillator; 



H 



f 



Ap 2 uJoAq 2 



+ 



2 2 

denotes the contribution from quantum fluctuations and 



(24) 



(25) 



(26) 



Hp = exp 



Aq 2 fd_" 
2 Uo 



(27) 



denots the contribution from the external perturbative potential. The external potential 
can be even or odd. Without loss of generality, we denote it as 



From Eqs. 



V{q) = KQ(q), 
and (|16|), we have 



9(g) 



sin(/c g), for odd 
cos(k q) 1 for even. 



(28) 



q = p, 

V = -^Iq - K eff e'(q)5 T , 
G = AUG, 



(29) 



n 



i 



8G 2 



-2U 2 



where 



K, 



eff 



K exp 



Tik 2 G^ 



(30) 



is called the effective potential, whose physical meaning will be discussed below. G' is the 
first derivative of with respect to q. It is clear to see from Eq. (p9|) that, in the time 
interval nT <t < (n + 1)T, the harmonic oscillator takes the free motion governed by 
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(31) 

col 
2 ' 

while, at the time t — (n + 1)T, it is kicked by the external potential. Therefore, before and 
after the kick the momentum and its fluctuation undergo jumps: 

Pn+1 (T + ) = p n (T-) - Ke-^^e'iqniT-)), 

n n+1 (T + ) = n„(T-) + -jKe-^^Q(q n (T-)). (32) 

Eqs. ( |3~TD and (|32|) are coupled differential equations. From these equations, we can readily 
see that the squeezed state dynamics differs from the classical one in two ways. (1) There 
are two additional dimensions for the squeezed state motion, namely, the particle moves in 
a four dimensional extended phase space. (2) The classical quantities (q,p) are coupled with 
the quantum fluctuations, which make the semiquantal motion complicated. On the one 
hand, because of these two additional dimensions in phase space, we expect that invariant 
curves in classical phase space would not be able to prevent the trajectory from penetrating 
or crossing them semiquantally. On the other hand, since the quantum fluctuations are 
always positive, they equal zero only in the limiting case of % = 0; the effective potential 
strength K e ff is always less than K. The reduction of the effective potential acting on 
the wave packet leads to the suppression of chaos. These two mechanisms coexist in the 
semiquantal system. They compete with each other and determine the dynamical behavior 
of the underlying system. Therefore, we expect that quantum fluctuations may not only 
enhance chaos but also suppress classical diffusion as well. This argument will be nicely 
illustrated in the following. 

In this section we restrict our calculations on the potential 

V(q) = -Ksm(q). (33) 

However, we should point out that the main conclusions given in this section do not depend 
either on the parity of the potential or the sign of K. 



P 
G 

ri 



-u q, 
4UG, 
1 

8G2 " 



2IT 



A. Enhancement of chaos 

In solving Eqs. fl3~T|) and (0), we used the seventh and eighth order Runge-Kutta formula 



with adaptive stepsize control. The permissible error is fixed at 10~ 12 . In Fig. 1(a), we plot 
the classical phase space (q c i,Pd) for a trajectory starting from (0,0) and evolving 10 4 kicks. 
The parameter K = 0.8 and a = where 

a = £> (34) 

is the ratio between the angular frequency of the kicks 27r/T; T is the period 

of the kicks) and the angular frequency of the harmonic oscillator u. In our calculations, 
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we put ujq = 1. It is obvious that, in classical phase space, regular (stable islands) and 
chaotic regions coexists. Fig. 1(b) shows the time evolution of expectation values (q,p) of 
the wave packet for 10 4 kicks. The wave packet starts from (q ,p ,Go,U ) = (0,0,0.5,0), 
with h — 0.1. The selection of initial conditions G = 0.5, and n = follows the minimum 
uncertainty and least quantum effect conditions given in Sec. II. Therefore, the initial wave 
packet has the same width in both q and p directions, that is, it is a coherent state. If 
there are no kicks, the wave packet starting from this point will evolve exactly along the 
classical particle's trajectories forever. The fluctuations both in momentum and coordinate 
keep constant, and are independent of time. In this case, the squeezed state dynamics 
exactly describes the classical one. Now, if we switch on the kick, the situation becomes 
quite different. As is shown in Fig. 1(a), the initial point just lies in stochastic seas, thus 
it is evident that in the classical case the trajectory will never enter into stable islands due 
to the existence of invariant curves. However, as we predicted the invariant curves are not 
able to prevent the trajectory from crossing it via other dimensions semiquantally. This is 
demonstrated by Fig. 1(b), where all stable islands in the classical phase space are "visited" 
by the semiquantal trajectory. 

As a quantitative verification, we have numerically calculated the maximal Lyapunov 
exponent 

A = lim A n (35) 

n— »oo 

for the trajectories in both cases. The time behavior of \ n is shown in Fig. 2. It demonstrates 
the coexistence of enhancement and suppression. At the initial stage the enhancement 
mechanism is dominant. However, after a certain time, A n of the squeezed state becomes 
larger than its classical counterpart, which means that the enhancement mechanism becomes 
dominant, and consequently leads to enhancement of chaos. Furthermore, the chaotic motion 
in the extended phase space is characterized by two positive Lyapunov exponents in four- 
dimensional phase space (q, p, hG, U) which could be verified readily. 

B. Suppression of chaos 

We would like go to another limit, namely, very large perturbation, to investigate the 
suppression of chaos. Classically, when K increases, the motion becomes more and more 
chaotic. For a sufficient large K such as K = 6, the classical motion is completely chaotic, 
as shown in Fig. 3(a), where a — Like Fig. 1, Fig. 3 is for a trajectory starting from 
the origin and evolving 10 4 kicks. The classical chaotic and diffusive process is easily seen 
from the evolution of this phase plot. To demonstrate the suppression of chaos (or diffusion 
process), we start a wave packet from (0,0,0.5,0) in the four- dimensional (4D) squeezed 
state phase space. The evolution is shown in Fig. 3(b). Comparing Figs. 3(a) and 3(b), 
it is obvious that, in the classical case, the phase space is chaotic and diffusive, whereas in 
the semiquantal case the diffusion process is largely slowed down and suppressed. There are 
invariant-curve-like structures that appear in the semiquantal phase space. These structures 
seem to form a barrier for diffusion and thus suppress chaos. The suppression of chaos is 
quantitatively demonstrated by the large decrease of A„, as is shown in Fig. 4, where the 
suppression mechanism is most important. 
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To illustrate the suppression, we plot variation of K e ff with time (in units of kicks) in 
Fig. 5. This plot indeed demonstrates that the effective perturbation strength is much 
less than its classical counterpart for most of the time during the evolution. This is the 
reason for the suppression. In fact, the deduction of the effective potential acting on the 
wave packet has a clear physical picture. The width of a wave packet centered at (q,p) 
in coordinate space is Aq = and the external potential has a wavelength of 2-rr/ko. 

Therefore, there are, in fact, m(— y/hGko/2ic) periods of external potential acting on the 
wave packet simultaneously. This is quite different from the classical model, where only one 
kick acts on the particle at one time. Since the external potential is negative in some places 
and positive in other places, the wider the wave packet, the larger the number of m, and 
therefore, the smaller the effective potential acting on the harmonic oscillator. However, if 
the wave packet Aq is so small that it is smaller than the period of the external potential, 
then the effective potential is large. As a matter of fact, the effective potential K e ff in 
Eq.(pO|) can be written as 

"2^J- (36) 

As a significant evidence of the suppression, it is convenient to calculate energy diffusion 
with time n (in units of kicks) for an ensemble of trajectories. The diffusion is defined by 
(E n ) subtracts initial averaging energy (Eq), where (• ■ ■) means ensemble the average over 
many trajectories. In our calculations we have taken such an ensemble averaging over 10 4 
initial points which are uniformly distributed inside a disk area centered at the origin of the 
phase space. For the classical one, E n = |(p^ + q^) c u an d for semiquantal dynamics, E n is 
defined by 

= \(jPl + + \n (^T + ^iGn + ^n) (37) 

In Fig. 6, we show the energy diffusion of K = 6 and o = 1/tt for classical and semiquantal 
cases. The suppression of classical diffusion is very obvious. 

C. Transition from enhancement to suppression 

We have seen so far that enhancement may happen at the small K regime, and suppres- 
sion at the large K regime. Now we would like to discuss the transition from enhancement 
to suppression by changing the strength of the external potential for fixed quantum fluc- 
tuations. Here we want to show that there exists a threshold value of K c distinguishing 
enhancement from suppression. 

To this end, we need to take an appropriate ensemble average over many trajectories 
in phase space. However, since the classical phase space of the KHO model is unbounded, 
it is impossible to do such an average over the whole phase space. This makes numerical 
works very difficult. After many numerical experiments, we find a compromise, namely, we 
take the average over a disk centered at origin with radius it. We spread 15 x 15 initial 
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points uniformly distributed inside this area. In classical case, we calculate the Lyapunov 
exponent for each trajectory after 10 4 kicks, and plot the averaged value denoted as (A^) 
in Fig. 7. This averaged value is in analogy to the Kolmogorov entropy in a bounded 
system. However, strictly speaking, this quantity cannot be called Kolmogorov entropy. 
Nevertheless, this parameter captures more or less chaoticity of the underlying system. In 
the case of semiquantal dynamics, since we have 4D extended phase space, we always have 
two positive Lyapunov exponents. We add these two values, and denote the result as (X sq ) ■ 
It is plotted in Fig. 7 in comparison with the classical result. 

From Fig. 7 we can draw the following conclusions: (1) There exists a certain threshold 
value of K c . Before this point, (X sg ) > (Xd), which means that the degree of chaos is 
enhanced; after this point, (X sq ) < (X c i), chaos is suppressed. This critical value K c changes 
with h. (2) At the region of K 3> K c , (X sq ) fluctuates around a certain value. It does not 
change with K. (3) The enhancement and suppression depends largely on Ti. 

The results discussed in this section are restricted to an irrational frequency ratio. One 
might ask, whether our conclusion also applies to the rational frequency ratio. It is well 
known that the KHO model is a degenerate system out of the KAM theorem. In classical 
phase space, there exists a slow diffusion along the stochastic web for any small value of 
perturbation. Our numerical results also show enhancement and suppression. We give one 
example of rational frequency ratio a = 1/4 and K = 6 in Fig. 8 for suppression. The 
corresponding Lyapunov exponent is shown in Fig. 9. 

Finally, we would like to say a few words about the initial conditions and parity of the 
external potential. We have performed a wide range of numerical investigations, and found 
that the above discussed qualitative and quantitative conclusions are independent of the 
selection of the initial condition and the parity of the external potential. However, the 
selection of the initial condition must be physically meaningful, as we discussed in Sec. II. 

Before concluding this section, we would like to discuss the connection of suppression to 
the dynamical localization. In fact, this is a challenge to the squeezed state approach for 
this subtle phenomenon. We argued that the dynamical localization observed in the kicked 
rotator is due to a suppression of the chaos discussed above. In fact, in the limiting case of 
ujq = 0, the KHO model Eq. [^] is reduced to the kicked rotator model, in which chaotic 
diffusion is completely suppressed by the quantum fluctuations and results in dynamical 
localization, a well established fact observed numerically by Casati et al almost 20 years 
ago, and confirmed recently by experiment ||17|| . This was nicely illustrated by Zhang and 
Lee H with the squeezed state approach. 



IV. DIFFUSION AND LOCALIZATION 

In Sec. Ill we showed how quantum fluctuations enhance and suppress chaos. This 
fact will definitely affect diffusion behavior. For instance, in the limiting case, when the 
suppression becomes dominant localization is expected to happen. In this section we shall 
give a detailed study of this. In particular, we concentrate on the large K regime. This 
is the most difficult region in pure quantum computation. As we shall see the squeezed 
state approach not only provides an easy way to do numerical calculations but also makes 
it possible to do some analytical estimations. 
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The energy E n of the kicked harmonic oscillator in the squeezed state approximation 
[Eq. (|37| can be written as two parts, 

E n = E c n + El (38) 

E^ contains the first two terms in Eq. (p7|), which is due to the motion of the centroid of a 
wave packet. They mimic the effect of classical diffusion (ECD). E? includes the last three 
terms in Eq. (|3~7D, and is attributed to the effects of quantum fluctuations (EQF). These 
two kinds of effects are the main ingredients of the diffusion process in the squeezed state 
dynamics. 

The ratio a is an important quantity, as we shall see soon. We take a, the golden 
mean value a g = (v^5 — l)/2, and its continued-fraction expansion r/s: 2/3,3/5,5/8, ■•• 
as examples, r and s are generated by the Fibonacci sequence defined by F — l,iq = 1, 
and F n = F n _ 2 + F n _i for n > 1. Without loss of generality, in all calculations, we keep 
parameters % = 1, K = 6, ko = 1, and uq = 1, and the initial point is chosen as (0, 0, 0.5, 0), 
which corresponds to the ground state of unperturbed quantum harmonic oscillator. 



A. Numerical results 

Figures 10 and 11 and 12 and 13 show our numerical results of energy diffusion of 

n 9 ) = {£ co t1 (39) 

{ Ksm(k q) 

for even and odd parity, respectively. Now we discuss these two cases separately. 

Even potential: In this case, for all rational frequencies, the energy will finally goes 
quadratically with time i.e. E n ~ n 2 . As shown in Fig. 10, the slope equals 2 asymptotically 
in the double logarithmic plot. However, for the case a = r/s with relatively larger r and 
s, the diffusion starts only after a certain time. Before this time the energy diffusion is 
localized. The transient time depends on the frequency ratios, and is approximately of the 
order of (er — (T g ) ■ We call this transient region a transient dynamical localization region. 

For the irrational case, dynamical localization occurs, as clearly demonstrated in Fig. 
11. This significant phenomenon has been observed and investigated in various quantum 
systems in past few years. 

It is worth pointing out that for two trivial cases, i.e. a = 1/1 and 1/2, our squeezed 
state results given here agree completely with the quantum analytical results of Ref. |2^ 



which has been the only existing analytical results of the quantum diffusion of this model 
up to now. This demonstrates the usefulness of the squeezed state approach. Moreover, 
by using the squeezed state approach, we have also recovered the quantum results obtained 



numerically by Borgonovi and Rebuzzini ||21|| . For more details, see Sec V. 



Odd potential: In this case, quadratic law is observed only in the case of rational frequency 
ratios a = r/s with odd s. For other situations the energy diffuses linearly with time 
approximately, see Figs. 12 and 13. 
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B. Analytical estimates 



The above numerical results can be understood by analyzing evolution equations 
In fact, we can analytically derive the energy diffusion by studying system fl2~9"|). Starting 
from Eqs. (fffi) , we find that when the external perturbation is absent, the two degrees 
of freedom (DOF's) (p, q) and (G, II) are decoupled, and each undergoes free motion. In 
terms of action-angle variable, the Hamiltonian of the free motions can be expressed as Eq. 
( Pjjp . From this formula, we have already seen that both free motions of the two DOF's 
are degenerate. It is this degeneracy which makes resonance between the two frequencies 
possible in phase space. Consequently, the squeezed state dynamical behavior of the kicked 
harmonic oscillator is quite different from that of the kicked rotator ||. That is, the motions 
of the centroid and the fluctuations of the wave packet behave like an oscillator with fixed 
frequencies ujq and 2ujq, respectively. 

However, when kicks are added, the two degree of freedom becomes coupled, and energy 
may start to diffuse. It is convenient to express the evolution of system in terms of action- 
angle variables. From Eqs. (|A4|) and (|A5|) one can readily obtain 4D maps 



4+i = h ~ Keff&' k J2I n /u; sm(j) n k Q J2I n /u} cos(j) n 



<Pn+i = <p n + u T + K eff e' \k y2I n /u sm(j) n ) k / y2I n u sin <f> n , 

J n+1 = J n + K eff jkl\l (4J n + l) 2 - 1 sin 9 n e (k yj2I n /u sin 0„) , (40) 

d n +i = d n + 2u T - K eff k 2 h ( 1 / AJn + ^= cos# n J 6 (k ^2I n /uj sin (j) n \ . 

(4</ n + l) 2 1 J 



With this 4D map, we are able to performe analytical estimate of the energy diffusion. We 
shall treat it at two different limiting cases. 

Classical diffusion effect (h = 0) 

In the classical limit case, % = 0, and E[ = 0, and the effect of classical diffusion becomes 
dominant. Therefore, the change of energy during one kick is 

AE c n = k KQ\k q n ) {p n cos u T - q n u smu T) + h 2 K 2 Q' 2 {k q n ). (41) 

For K 3> 1, the orbit can be supposed to be approximately ergodic. After ensemble averaging 
over variables p and q, the first two terms vanish approximately, thus we obtain linear energy 
diffusion, 

El ~ (AE c n )n « \k 2 K 2 n. (42) 

Note that the average of the first two terms of A-E£, though is much smaller than the last 
term for large K, is nevertheless not exactly zero. This results in some oscillations of E n 
around the linearity (see Figs. 12 and 13). 
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Effect of Quantum Fluctuations 



For the second limit case, suppose a wave packet starts from (q,p, G, II) = (0, 0, 0.5, 0) with 
an even potential, and the center of the wave packet keeps fixed, thus E n = 0. The energy 
diffusion is caused purely by the effects of the quantum fluctuations. In this case we shall 
analyze the diffusion process for two different frequency ratios, i.e., rational and irrational. 
For the rational frequency ratio, let us take the simple case of a = 2/3 as an example. During 
a time of 3T, there are three kicks acting on the harmonic oscillator. Since the frequency of 
fluctuation is 2u , (G, II) evolves four periods. Note that the effective amplitude of a kick 
acting on the wave packet is K e ff rather than K. Among these three kicks, only that one at 
a relative small G affect the free motion of the oscillator significantly. We call this kick the 
effective kick, the effects from other two kicks can be neglected due to a very large G, and 
consequently a very small K e ff. At the time when the next effective kick is in action, G is 
approximately the same because of the resonance, see Fig. 14. Therefore, the increment of 
II is almost constant, which means that n 3ra »s n. Thus, from Eq. Q3"7|), we obtain 

E n « n 2 , (43) 

which gives rise to the quadratic law observed in Figs. 10 and 12. 

If a is an irrational number, a very interesting thing will happen. From Eq. ( |A3| ) we 
know that the angular variable when a kick is added is 

9 n = 2nna + o (mod2ir). (44) 

This is nothing but a pseudorandom number generator, indicating that the jump of II may 
happen in upper (II > 0) and lower parts (II < 0) with the same probability, and thus 
the increment of energy in the upper part will be canceled out by the decrease in the lower 
part. This leads to the localization phenomenon observed in Fig. 11. This localization 
mechanism, resorting to a pseudorandom number generator, reminds us what happens in 
the kicked rotator, where the localization is related to Anderson's localization for a quantum 
particle propagating in a one- dimensional lattice in the presence of a static-random potential 
[ill . Our results imply that it might also be possible to construct a connection between the 



kicked harmonic oscillator and the Anderson's problem in the framework of the squeezed 
state approximation. The mechanism discussed can also explain the transient dynamical 
localization that occurs in the case of the rational frequency ratio, as shown in Fig. 10. 
Since during the time t < (er — cr s ) _1 , a rational number behaves just like a pseudoirrational 
number, a transient localization phenomenon occurs. 



General case 



As to the general case of system (|29|) , both the ECD and EQF may coexist. To illustrate 
this, we consider the case of a = r/s, where r and s are coprimed integers. Suppose s is 
odd. As we have explained above, between two effective kicks, (q,p) and (G, II) evolve freely. 
Thus the angle variables of (q,p) at the two successive effective kicks are and <fi + 2wr, 
respectively, and that of (G, II) are 6 and 6 + Anr. From Eqs. ( |A5|) and (p9|), we find that in 
this case, the increment of Ap and All have the same sign, which means that both the ECD 
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and EQF are excited, which is independent of the potential parity. Because the diffusion due 
to the EQF is ~ t 2 , which is much faster than that of the ECD (~ t), thus asymptotically 
t 2 diffusion shows up as shown in Fig. 12. However, if s is even, there is one additional 
effective kick between the above mentioned two, namely, the s/2th kick, at which the angle 
variable of (G, II) is 9 + 2nr, while that of (q,p) is + nr. Thus, for the even potential 
case, the changes of II due to the two consecutive effective kicks have the same sign, which 
implies that the EQF is excited and t 2 diffusion will show up. For the odd potential case, 
the change of II due to the two consecutive effective kicks have opposite sign, the EQF is 
thus suppressed, and we obtain the linear diffusion as seen in Fig. 13. 

For the irrational a case, the EQF is suppressed and the ECD becomes dominant. If 
the (p, q) happens to be a fixed point in (p, q) plane as is the case of Fig. 11, localization 
occurs. This is the reason for the different diffusion behavior of irrational o in Figs. 11 and 
13 for even and odd external potentials, respectively. Please note that (p, q) = (0, 0) are the 
expectation values of all the eigenstates of the harmonic oscillator, thus the localization we 
observed in Fig. 11 is not restricted to the case of the ground state, as we discussed up to 
now, but is very general. 

C. Transition from localization to diffusion 

The results discussed above focused on diffusion at very large perturbation, in this case 
the underlying classical system is completely chaotic. As a further example we would like 
to demonstrate a very interesting and important phenomenon in quantum mechanics, i.e. 
the tunneling effect. We start a wave packet from point (0,7.5) in the classical phase space 
(q,p). Here we have a = 1/5 and V(q) = Kcosq with K = 0.5. The wave packet has 
parameter Go = 1/2 and n = according to the minimal uncertainty principle. It is a 
Gaussian wave packet. The classical phase space is shown in Fig. 15(a). We see that the 
starting point lies inside a stable island. Classically, a trajectory that starts from this point 
will never be able to escape. However, in the quantum case, the situation becomes very 
different. We expect that if the width of the wave packet is much smaller than the size of 
this stable island, the wave packet will be confined by this stable island, and thus lead to 
localization. However, if the wave packet becomes wider than the size of the stable island, 
it will spread out. Here we demonstrate this quantum phenomenon by the squeezed state 
approach. In Figs. 15(b)-15(d) we plot the energy evolution for different Planck constant 
ft, which corresponds to different widths of the initial wave packet. At ft = 1 and 2 we 
observed a localization phenomenon as in other quantum systems. The energy oscillates 
around a certain value. When we increase ft further to ft = 5. A transition from localization 
to derealization occurs which is shown in Fig. 15(d). 

V. COMPARISON WITH QUANTUM RESULTS 

To give the reader a clear picture about the accuracy of the squeezed state approach, 
we would like to compare our results with those obtained from pure quantum computation. 
However, as mentioned above, since diffusion occurs in the whole unbounded phase space 
and cannot be reduced to motion on a cylinder like the case of the kicked rotator, a pure 
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quantum (numerical) investigation is very difficult, in particular, in the large K regime. 
Nevertheless, with a large amount of CPU time, one would be able to obtain some results 
in the small K regime for short time evolution. This is why there are only limited available 
quantum results in this regime. We will compare them with our squeezed state results here. 

First, let us look at the results obtained in Ref. f20fl , already shown in Fig. 10. As for 
the case of 1/1 and 1/2, our squeezed state results agree completely with the quantum one 
in Ref. [2T| . In this case both squeezed state and pure quantal analyses predict t 2 growth 
of the energy. We should mention that it was conjectured in Ref. p0[ that for the general 
case of 1/q, the energy growth should less than t 2 . From our squeezed state analysis we 
concluded that there are only three different diffusions: t 2 , t and localization. Therefore, 
our squeezed state analysis also agrees with Ref. pOj's prediction. 



Now we turn to the results obtained by Borgonovi and Rebuzzini pTJ . Since the time unit 



given in their pictures is not clear, we are not able to make any quantitative comparison 
with our squeezed state results. Therefore, we performed quantum calculations by using 
our own program (see Appendix B). All the system parameters are kept the same as that 
used by Borgonovi and Rebuzzini. The results are given in Figs. 16 and 17. These pictures 
correspond to different diffusion behaviors. In Fig. 16, at 2 diffusion is obtained (t is in units 
of kicks); the squeezed state approach also gives rise to a t 2 diffusive behavior, although there 
is a difference in prefactor. In Fig. 17, both the quantum and squeezed state results show 
localization around approximately the same energy. Please compare these two pictures with 
Figs. 1 and 8 of Borgonovi and Rebuzzini f2"lfl , respectively. 

Our quantum computation technique (Appendix B) is different from that used in Ref. 
[2(J. For a self-consistent test, we have used the same parameter as that of Fig. 2 in Ref. 
1 2(1, and computed the energy diffusion with our method. We found that our results agree 



with those of Ref. |2(J in every detail. 

As already emphasized, because of the unbounded phase space, the quantum computa- 
tion is very time consuming, even for small perturbation. For instance, about ten days CPU 
time (IBM RISC System/6000 42T, with 192 Mbyte RAM) has been spent for Fig. 16, and 
20 days CPU time for Fig. 17. 



VI. CONCLUSIONS AND DISCUSSIONS 

Applying the squeezed state approach to the kicked quantum harmonic oscillator, we 
illustrate how the quantum fluctuations affect the classical dynamics. We have shown that 
chaoticity can be enhanced as well as suppressed by the quantum fluctuations. A transition 
from enhancement to suppression is observed when we change the strength of the kicks. 

Moreover, with this squeezed state approach, we are able to investigate the energy diffu- 
sion. Three different energy diffusion have been observed for the kicked quantum harmonic 
oscillator, namely, localization, linear diffusion, and quadratic diffusion. The localization is 
due to strong suppression. 

Though it is a kind of approximation, the squeezed state can mimic many true quantum 
behaviors such as those demonstrated in this paper. Moreover, there are examples showing 
that the squeezed state approach gives rise to better results than the semiclassical method. 
However, under which condition or how far the squeezed state approach can go beyond the 
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semiclassical method is still an open problem that deserves further numerical as well as 
theoretical study. 

Finally, we would like to remark on the quantization of a quantum system whose classical 
counterpart is chaotic. As is well known, this is a tough problem which has attracted tremen- 
dous attentions in last two decades. Among many others, Gutzwiller's trace formula might 
be the most plausible one |p7| . However, this approach encountered difficulty of divergence, 



although many important contributions have been made to overcome this difficulty. We 
believe that the squeezed state approach might be an alternative way that can contribute 
to this. Recent successful application of this method by Pattanayak and Schieve [FL\ to 
calculate eigenenergies of a chaotic system sheds light on this direction. 
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APPENDIX A: QUANTIZATION OF THE HARMONIC OSCILLATOR BY THE 

SQUEEZED STATE APPROACH 

In this appendix, we would like to demonstrate how the squeezed state works when 
applied to a simple quantum system, a harmonic oscillator, which is the KHO system with 
a zero external potential. The harmonic oscillator is a simple but very important model in 
quantum mechanics. It is an integrable system, and its eigenenergies as well as eigenf unctions 
can be obtained analytically. Therefore, this model is very suitable for testing approximate 
methods such as the WKB method and others. 

It is well known that the WKB approximation can give us exact eigenenergies for this 
integrable system. However, it cannot yields exact eigenfunctions, in particular, for the 
low-lying eigenstates. It gives only the envelope of the wave function in the semiclassical 
limit h — > 0. In this appendix, we shall demonstrate that when applying the squeezed state 
approach to the harmonic oscillator model, one can obtain the energy levels precisely as well 
as the eigenfunctions. 

The harmonic oscillator has the Hamiltonian operator 

* = ? + 1r < A1 > 

Applying the squeezed state approach to this system, one can easily obtain 

H = ^ + ^f + l(^G + + • (A2) 
This Hamiltonian can be expressed in terms of action-angle variables 

H = u Q I + 2u {J+ti~), (A3) 



17 



where 



I=^<fpdq, J=±.fud(hG). (A4) 
The transformations between (q,p), (G, IT) and (1,0), (J, 9) have the following forms: 



21 /u Q sin< 



p = y 21oj cos ( 
1 /, . 1. 



G = — l(2J + -) - J2J(2J + l)cos6) , (A5) 
too V 2 v 



f ^2J(2J + l)sin0 
(2J+ |) - ^/2J(2J + 1) cos 0* 

From Eq. (|A5|) one can see that the motion of both degrees of freedom are degenerate, namely, 
^f- is independent of / and ^y- independent of J. Furthermore, (G, II) are decoupled from 
(q,p). Thus the centroid of the wave packet goes exactly along the classical trajectory. 
While the fluctuations in momentum and position are time independent, i.e., the width 
of the wave packet keeps constant. From the minimal uncertainty principle for the initial 
condition mentioned in Sec. II, we have G = ^ and II = 0. Therefore, the wave packet 
along the periodic orbit always keeps its form as a coherent state. 

The time evolution of both (q,p) and (G, IT) are periodic with period T (= 2tt/cjo) and 
To/2, respectively. So we can apply the EBK quantization to the extended phase space, 

I = nh, J = mh, n, m — 0, 1, 2, • • • , (A6) 

Substituting I and J into Eq. (|A3|) and keeping in mind that G = ^ and II = 0, and thus 
m=0, we obtain the squeezed state eigenenergy, 

E n = hu {n + h. (A7) 

This is exactly the eigenenergy of the harmonic oscillator. Here the zero point energy ^Huq 
comes into the formula in a very natural and straightforward way. This is quite different 
from the WKB method. In the WKB method, comes from the Maslov phase correction 
which is necessary because of the singularity of the wave function. In the squeezed state 
approach, since we do not have any singularities, the Maslov-Morse correction is incorporated 
by the extended variables G and II. Furthermore, the energy of the system is in the form of 
the expectation value of the underlying Hamiltonian operator, whereas in the usual WKB 
method, the energy is taken to be the classical form, i.e., E = H c \. This is one of the reasons 
that people call this approach the semiquantum approach. 

Let us now construct the eigenfunction by the squeezed state approach. We see that 
when the trial wave function is transformed to 

!*(*)> = exp (^) |¥(f)>, (A8) 
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the derived variational equations of motion remain invariant. Substituting \^{t)) into the 
Schrodinger equation, we can determine the equation of X(t), 

X(t) = jT dt'^{t')\ih^ - H\9tf)) = \ G + A D . (A9) 

The second part of the integral corresponds to the dynamical phase. We denote it Ad. The 
first term is the geometrical phase noted as Xq, which is 

\ G = - [ (pq - qp)dt + H f tlGdt. (A10) 
2 Jo Jo 

Since the motions of (p, q) and (G, II) are periodic, this geometrical phase is the Aharanov- 
Anandan form of Berry's phase. During the evolution, each points along the periodic orbit 
acquires a phase factor. However, the dynamical phase does not change during the evolution, 
only the geometrical phase matters. So the eigenfunction is a weighted sum over points of 



the commensurate periodic orbit flO|| . The weight factor at each point is an appropriate 
geometrical phase. Furthermore, as mentioned above according to the requirement of the 
initial condition, the initial wave packet is a coherent wave packet, and it does not change 
its form when cycling along the periodic orbit. Substituting expressions for p and q in Eq. 
into Eq. ( |A10|) , and keeping in mind that II = 0, we can easily evaluate the integral 



and obtain the geometrical phase at time t, which is 

X G (t)=nh(p, (All) 

where = 7r/2 — rj. Thus the eigenfunction for the bound state having the eigenenergy E n 
in Eq. ([AT]) is 

c e l -^\a)d(f) = C e- m i\a)dr], (A12) 
Jo Jo 

where \a) is the coherent state, is the angle in the (uioq, p) plane, and C is the normalization 
constant. This is nothing but the number eigenstate \n) given in Eq. (0) except for the 
prefactor. This constant C can be easily calculated by the normalization. Therefore, in the 
coordinate representation, the wave function is 

^n(q) = C r d V e- m *>(q\a) (A13) 
Jo 

This is the exact wave function of the harmonic oscillator. 



APPENDIX B: PROCEDURE OF QUANTUM COMPUTATION 

In this appendix, we describe our procedure of quantum computation. Since the Hamil- 
tonian is periodic in time, the Floquet theory can be applied. The time evolution can be 
reduced to evolution of the eigenstate over one driving period, 

Mt + T)) = U(T)\*(t)), (Bl) 
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where 



U(T) = UfreAick = exp (-i^p) ex P ^ B2 ) 

is the Floquet operator. 

To simulate quantum diffusion in this degenerate system ( p2j ), the Fourier spectral 
method is employed. The time interval of free propagation is divided into many slices, 
each having a width of A. For each slice the evolution operator is factored into a product of 
kinetic and potential propagator arranged in a symmetric way so that a full potential step 
is sandwiched between two half kinetic steps, namely, 

- (-*¥) - - M£) - (-^) - ) + °' a3 » <-> 

This technique of symmetrically splitting the kinetic propagator reduces the error introduced 
by neglecting the commutator between the kinetic and potential operators. The error is 
reduced to 0(A 3 ) from 0(A 2 ) in a nonsymmetric splitting. The kinetic propagation is 
carried out in momentum space, since in this space the time evolution is simplified as 
multiplication. The potential step is performed in coordinate space for the same reason. 
The kick step, performed once per period, is also done in coordinate space. A fast Fourier 
transform routine is used to transform wave function between these two spaces. 

Since the KHO model is a degenerate system, a wave packet may diffuse rapidly, even 
to infinity in both coordinate and momentum space. The average energy of the wave packet 
may reach a rather high value during the diffusion. This amplifies the error caused by the 
approximation made in Eq. (p3|). Therefore, the self-adaptative procedure is used to adjust 
the time slice in Eq. ( |B3[ ) in each period to make sure that the width of the time slice 
is much smaller than the inverse energy. Second, both coordinate and momentum spaces 
should be large enough. So a large number (32 768) of Fourier components are used in our 
computations. 
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FIGURES 



FIG. 1. Classical phase space (a) and time evolution of expectation values (p, q) of a wave 
packet (b) at K = 0.8 for an irrational frequency ratio a = 1/tt. One classical trajectory 
starts from (0,0). The wave packet starts from a point having an initial fluctuation parameter 
(Go,IIo) = (0.5,0), and h = 0.1. The initial wave packet is a coherent state, i.e. it has same width 
in both p and q direction. 

FIG. 2. Time (in units of the kicks) behavior of \ n for the trajectory shown in Fig. 1. The 
increment of \ n after a certain time in the semiquantal case indicates the enhancement of chaos. 

FIG. 3. Same as Fig. 1 but for K = 6 with an irrational frequency ratio a = 1/tt for 
classical (a) and squeezed state (b) cases (h = 1). The semiquantal phase space shows an obvious 
suppression of the classical diffusion. 

FIG. 4. Time behavior of X n for the trajectory shown in Fig. 3. The large decrement of A n in 
the semiquantal case demonstrates the strong suppression of classical diffusion. 

FIG. 5. Time evolution of the effective external potential K e ff for the orbit shown in Fig. 3. 

FIG. 6. Energy diffusion at K = 6 for an irrational frequency ratio a = 1/tt for classical and 
semiquantal cases. The ensemble averaging is taken over 10 4 initial points which are uniformly 
distributed in an area of the disk centered at (0,0) with a radius of ir in the phase space. The 
diffusion coefficient in the semiquantal case is obviously much smaller than that of the classical 
case, which indicates a strong suppression of classical chaos. 

FIG. 7. Transition from enhancement to suppression. The averaged Lyanpunov exponent vs 
the external potential K for classical and semiquantal cases with h = 0.1 and 1. The average is 
taken over 15 x 15 points uniformly distributed inside a disk of radius tt centered at (0,0). Here 
a = 1/ir. 

FIG. 8. Demonstration of suppression for a rational a = 1/4 frequency ratio and K = 6. (a) 
Classical phase space, (b) Semiquantal (q sq ,p S q) with U = 1. 

FIG. 9. The Lyapunov exponent for the trajectory shown in Fig. 8. 

FIG. 10. Evolution of the energy with time (in units of kicks) for the case of an even potential 
with a rational frequency ratio. Numbers in the plot indicate the frequency ratio. The case of 
1/1 coincides with that of 1/2. Note that the curves have slope 2 asymptotically, and a transient 
dynamical localization phenomenon shows up. 
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FIG. 11. Same as Fig. 10, but for an irrational frequency ratio a = — l)/2. The 
localization is obvious. 

FIG. 12. Same as Fig. 10, but for the case of odd potential. The dashed line with slope 1 is 
drawn to guide the eye. 

FIG. 13. Same as Fig. 11, but for the case of odd potential. The dashed line with slope 1 is 
drawn to guide the eye. 

FIG. 14. Evolution (G, U) plot for K = 6 and a = 2/3. 

FIG. 15. Transition from localization to delocalization of a wave packet driving by the quantum 
fluctuations. The wave packet starting from a stable island (qo,Po) = (0, 7.5). a = 1/5, cjq = 1 and 
K = 0.5. (a) Classical phase space, (b)-(d) Semiquantal energy diffusion, (b) h = 1. (c) U = 2. 

(d) n = 5. 

FIG. 16. Comparison of semiquantal (dashed line) and quantum (solid line) diffusion for 
a = 1/4, loq = 2/-7T, uj\ = 8tt, K = 0.5, and fi = 1. The trajectory starts from (qo,Po) = (3.15,0). 
The semiquantal trajectory has the same (qo,Po) an d (Go,ITo) = (^,0). Please compare it with 
Fig. 1 in Ref. [21]. 

FIG. 17. Comparison of semiquantal (top) and quantum (bottom) diffusion for a = (\/5— l)/2, 
loq = 1, K = 1, and h = 1. The trajectory starts from (qo,Po) = (15, 0). The semiquantal trajectory 
has the same (qo,Po) and (Go,IIo) = (^-,0). Please compare it with Fig. 8 in Ref. [21]. 
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